Full two-electron calculations of antiproton collisions with molecular hydrogen 
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Total cross sections for single ionization and excitation of molecular hydrogen by antiproton impact are 
presented over a wide range of impact energy from 1 keV to 6.5 MeV. A nonpertubative time-dependent close- 
coupling method is applied to fully treat the correlated dynamics of the electrons. Good agreement is obtained 
between the present calculations and experimental measurements of single-ionization cross sections at high 
energies, whereas some discrepancies with the experiment are found around the maximum. The importance 
of the molecular geometry and a full two-electron description is demonstrated. The present findings provide 
benchmark results which might be useful for the development of molecular models. 



PACS numbers: 34.50.Gb,25.43.+t 

A central point of atomic and molecular physics is the de- 
scription of charged particles moving in a Coulomb field. One 
of the simplest and most basic systems which provides an in- 
sight into dynamic processes of charged particles is the colli- 
sion of antiprotons with atoms. The heavy mass of the antipro- 
ton allows, first, for a semiclassical theoretical approach and, 
second, for the investigation of "slow" ionizing collisions. In 
contrast to positively charged projectiles, for antiprotons there 
is no complication from charge transfer. 

Further attention is drawn to this topic due to the upcoming 
FAIR 1 1] facility with the international collaborations FLAIR 
[0] and SPARC [3], both intending to investigate antiproton 
driven processes and even kinematically complete collision 
experiments. However, the design of FLAIR already requires 
a reliable knowledge of low-energy antiproton cross sections 
of residual gases as, e.g., molecular hydrogen. These experi- 
mental efforts complement the recent intensive studies on an- 
tihydrogen at CERN aiming to test the CPT invariance and to 
disclose the nature of antimatter gravity. 

Over the last decades a remarkable progress in the under- 
standing of interactions between antiprotons and atoms has 
been achieved (cf. and Refs. therein). The theoretical de- 
scription concentrated mainly on hydrogen and helium atoms. 
It was relatively easy to establish a full treatment of the for- 
mer target which only consists of one electron and one nu- 
cleus. For helium atoms, on the other hand, the dynamic 
electron-electron correlation effects turned out to be decisive; 
requiring much larger efforts for their correct description. Due 
to the above mentioned favorable properties, antiproton colli- 
sions on helium atoms became a benchmark system for study- 
ing electron correlation in atoms stimulating a large number 
of calculations which employed various theoretical methods. 
During the last ten years close-coupling calculations using ei- 
ther a spectral or spatial expansion of the two-electron wave- 
function 1110,0,0 yielded the most precise results where the 
latter usually take advantage of large-scale computing facili- 
ties. They have provided cross sections for single and double 
ionization which are mostly in agreement with experiment for 
intermediate to high impact energies while discrepancies still 
persist for low energies. 

Antiproton collisions with molecules have been studied ex- 
perimentally in a similar way as atoms concentrating mainly 



on ionization cross sections |0] and stopping powers JK 
where for the latter rather diverse results were obtained at 
low energies. In contrast, the theoretical work on collisions 
involving molecules is still comparably sparse (cf. Ill 111 and 
Refs. therein). Certainly, the description of a four-particle 
system like a hydrogen molecule, consisting of two electrons 
and two nuclei, is a further step in complexity compared to 
a helium atom. Recently, the ionization and excitation cross 
sections ifTTl fl2tl as well as the stopping power jl3ll for an- 
tiproton impact on molecular hydrogen were calculated us- 
ing spherical one-electron models for the hydrogen molecule 
IU4I1 . They could mostly reproduce the experimental antipro- 
ton results for impact energies E > 90 keV. The findings 
suggest, however, that for lower energies molecular as well 
as electron-electron correlation effects are important and have 
to be considered. An earlier work by Ermolaev [15] turned 
out to be unsatisfactory for E < 200 keV reproducing rather 
atomic than molecular hydrogen. Furthermore, two calcula- 
tions, both treating the target as a molecule, were performed 
by Sakimoto |[l6Tl and recently by the present authors lfl7ll for 
molecular hydrogen ions. It was shown that the calculation 
of only three orientations of the molecular axis with respect 
to the projectile trajectory are sufficient to obtain the ioniza- 
tion cross section 01711 . Currently new experimental data for 
antiproton collisions with molecular hydrogen are produced 
using the AD facility at CERN d . 

In response to the renewed experimental activity and the 
limited theoretical understanding a full two-electron close- 
coupling method has been developed. Converged cross sec- 
tions for single ionization and excitation of molecular hydro- 
gen are provided over a wide energy range from 1 keV to 6.5 
MeV on a dense energy grid. They demonstrate the impor- 
tance of a full two-electron description and of the molecu- 
lar geometry including different orientations of the molecular 
axis as well as the differences between atomic and molecular 
targets. To the best of the authors' knowledge no two-electron 
description for antiproton impacts on molecular targets has 
been introduced before in this energy range. 

The collision process is considered in a non-relativistic 
semi-classical way using the impact parameter method (cf. 
ll20ll ) which is known to be highly accurate for impact energies 
E > 1 keV [6]. The quantum-mechanically treated electrons 
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are exposed to the Coulomb potential of the molecular nuclei 
as well as the heavy projectile. The latter is assumed to move 
on a straight classical trajectory R(t) = b + vt given by the 
impact parameter b and its velocity v while t is the time. 

In the Born-Oppenheimer approximation the total wave 
function of the hydrogen molecule separates into the product 
(atomic units are used unless stated otherwise) 



r 2, Rnu) 



(k) 

where X \j are the eigenfunctions of the molecular vibra- In accordance with d, the transition probability 
tion, YF 1 the spherical harmonics, and 51 = (y, j, m) denotes 



-l} m (e,$)V fc (ri,r 2 ;i? nu ), 
(1) 



sum of products between the CI coefficients of ip k and ipj an d 
the one-electron matrix elements — calculated as in IU7I1 — be- 
tween the orbitals of the corresponding CI configurations. 

The coupled differential equations in Eq. © are integrated 
in a finite z-range —50 a.u. < z = vt < 50 a.u. with the 
initial conditions Cfc[R(to=— 50 /v)] — S k o, i.e., the target is 
initially in the electronic ground state ipQ with energy eo- The 
probability for a transition into the final state ipk at t/ = 50 /v 
for a fixed R nu is given by 



p k (b,E;R nn ,e^) = \c k (b,v,t f ;R nul Q,<S>)\ 



the vibrational and rotational quantum numbers. R nu = 
(i? nu , 0, $) and ri 2 are the position vectors of the nu- 
clei and the electrons, respectively. The wave function 
V'fc( r i) r 2! -Rnu) satisfies the electronic part of the time- 
independent Schrodinger equation (H e ) for an unperturbed 
molecule at a fixed internuclear distance R nu . The ip are ob- 
tained together with the eigenenergies e in a full configuration- 
interaction (CI) calculation l2ll I22I1 . The two-electron con- 
figurations are constructed from correctly anti-symmetrized 
products of orbitals which are eigenstates of the molecular 
hydrogen ion. The orbitals were obtained as in Ref. [17], 
where the radial part is expanded in B splines and the angu- 
lar part in spherical harmonics. More details on the extension 
of the close-coupling method from one-electron lfl7l fl8tl to 
two-electron targets are provided in 1I23I1 . 

For a fixed R nu the fully correlated wave function of the 
two-electron target molecule interacting with the antiproton is 
obtained by the evolution of the time-dependent Schrodinger 
equation in real time, 

i—^(r 1 ,r 2 ,t) = (ff e + M„t(ri,r 2 ,£)) *(n,r 3 ,*), (2) 

where the time-dependent interaction between the electrons 
and the projectile with charge Z p is expressed by 



Vint(ri,r 2 ,i) 



Z 



Z n 



| ri -R(*)| |r 2 -R(i)| 



(3) 



The interaction of the projectile with the nuclei, which leads 
only to an overall phase, is not considered here. 
The time-dependent scattering wave function 



#(ri,r 2 ,t) 



k 



(4) 



is expanded in the normalized eigenstates ip k of H e . Employ- 
ing this expansion in Eq. (0 and projecting with ipk leads to 
the usual set of coupled equations for every trajectory R(t), 
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Cj ( tpk I Vint I ipj } e 



(5) 



The two-electron interaction matrix elements in Eq. (0 can — 
according to Eq. © — be expressed as a sum of one-electron 
matrix elements in the orbital basis. The full two-electron in- 
teraction matrix element between tpk and ipj is therefore the 



Pk(b,E) 



(6) 



x^(iU)*7"(e,$)p (7) 

x Pk (b, E; R an , 8, $) sin edi? nu d9d$ . 



becomes orientation-independent by integration over R„ 
The corresponding cross section 



a k (E) = 2Tr J Pk {b,E)bdb, 



(8) 



can then be obtained by integration over b as is done for 
atomic targets which are spherical symmetric. The total cross 
sections for ionization <7i 0n and bound-state excitation cr cxc of 
the target can be obtained by summing up all partial cross sec- 
tions <7fc into states k [as given in Eq. ([8]l] with e k > Ik 2 and 
eq < e k < Ih 2 > respectively, where 7h 2 is me fi rst ionization 
potential of the hydrogen molecule. 

In this work, however, p k is approximated by 



Pk 



Pk(0,0)+p k -,0 + Pk 



7T 7T 

2' 2 



(9) 



using only the probabilities p k (0, <£>) for three orthogonal ori- 
entations of the internuclear axis with respect to the trajectory 
as is discussed in detail in Ref. 11711 . It was shown that this 
approximation is equivalent to Eq. (Q, if the integration is 
performed with a six-point quadrature formula Il24tl . Further- 
more, excellent agreement between the results using Eqs. (0 
and (0 were obtained for antiproton collisions with molecular 
hydrogen ions lfl6l [l7ll . 

The dependence of the ionization cross section on the inter- 
nuclear distance R nu in antiproton collisions with the molec- 
ula r hy drogen ion and the hydrogen molecule was examined 
in 01611 and 111 111 , respectively. In both cases an approximately 
linear dependence on R nu around the expectation value (R nu ) 
of the ground state was observed leading to Franck-Condon 
results which were found to be very close to the exact cross 
sections obtained by an integration over R nu . In this work, the 
calculations are accordingly performed for (R nu ) — 1.4478 
a.u. The two-electron basis consists of 3080 singlet states 
with total azimuthal quantum numbers M = 0, . . . , 4 and a 
maximum energy of e w 25 a.u. 

In Fig. [TJ the results of the time-dependent close-coupling 
calculations for (a) single ionization and (b) bound-state ex- 
citation of the hydrogen molecule by antiproton impact are 
presented. The present results are compared to the data avail- 
able in literature. The single-ionization cross section in the 
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Figure 1 : (Color online) Cross sections for (a) single ionization <Ti on 
and (b) excitation <t cxc by antiproton impact. Black solid curve with 
pluses, present results for molecular hydrogen; red solid curve with 
X, H2 model lllll ; blue dashed-dotted curve, scaled hydrogen atom 
H aC ai mil ; green doubly-dotted-dashed curve, hydrogen atom flUl ; 
maroon dotted curve, two times hydrogen atom ITHl ; violet dashed 
curve with triangles, H sca i by Ermolaev llal : green squares, CERN 
94 d ; red circles, CERN 90 (H . 



top panel is in excellent agreement with the experimental mea- 
surements for energies above 85 keV except for the data points 
at 500 keV. Below 85 keV the experimental data show a small 
discontinuous step and increase to a higher maximum than the 
present results which is situated in both cases around 40 keV. 
Note, in the extensive convergence studies performed in this 
work an enlargement of the basis always led to smaller values 
of the maximum. Below the maximum the experimental data 
fall off steeply in a similar way as the data for helium which 
were measured in the same occasion [91]. For helium, how- 
ever, the two lowest energy data points were withdrawn after 
a recent remeasurement ||4J]. The currently produced experi- 
mental hydrogen molecule data for low impact energies may 
help to clarify the trend below the maximum. 

The literature results obtained using a model potential and 
a hydrogen atom with scaled nuclear charge Z = 1.09 
il 1 1 411 are able to approximate the present calculations 
for energies above 50 and 100 keV, respectively. Though, they 
are throughout lower than the latter for these energies. For 
lower energies the models yield evidently too large cross sec- 
tions and in both cases show rather an atomic than a molecular 
slope by what they reveal their atomic nature. Below the max- 
imum also the lack of electron-electron correlation effects can 
be expected to become severe as in the case of the helium 
atom I4l 12311 . The calculations by Ermolaev 01511 . using also a 
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Figure 2: (Color online) Comparison of single-ionization cross sec- 
tions for antiproton impact on different targets. Black solid line with 
pluses, present results for molecular hydrogen; red dashed curve with 
circles, helium atom [ 23]; blue dot-dashed curve with stars, molecu- 
lar hydrogen ion 1 17J] (times 4); green doubly-dot-dashed curve with 
squares, hydrogen atom 1 18]. 



scaled hydrogen atom, are not satisfactory, since they follow 
for intermediate energies rather the data for a hydrogen atom 
multiplied by a factor 2. 

The lower panel of Fig. Q] compares the present excitation 
cross sections for molecular hydrogen with the existing liter- 
ature, i.e., the two already mentioned models lfTll[T8ll . Obvi- 
ously, the scaled hydrogen atom is not capable of reproducing 
the excitation cross section for molecular hydrogen despite its 
reasonable results for ionization for E > 100 keV. The model 
potential, on the other hand, is again an excellent approxi- 
mation for energies above 50 keV. This might have been ex- 
pected since the bound state energies and oscillator strengths 
of the model potential were found to be in good agreement 
with those of the hydrogen molecule in contrast to the ones 
predicted by the scaled hydrogen atom [14]. Note, for ener- 
gies above the maximum the cross section for excitation of the 
hydrogen molecule is quite similar to that of atomic hydrogen 
while for ionization it is rather comparable to twice the cross 
section of the hydrogen atom. 

In Fig. |2] the single-ionization cross section of molecular 
hydrogen by antiproton collision is compared to results for the 
helium [23] and hydrogen atom [18] and the molecular hydro- 
gen ion [17], where the latter is scaled by a factor 4. The com- 
parison shows that below the maximum the curve for molec- 
ular hydrogen decreases much faster with decreasing energy 
than is the case for the hydrogen and helium atom. The molec- 
ular hydrogen curve is on the other hand qualitatively similar 
to that of the molecular ion. That is, especially for these ener- 
gies the qualitative shape of the ionization cross section seems 
to be different for atomic and molecular targets. At low ener- 
gies ionization occurs mainly in a small region close to the 
nuclei where the electronic density and the expectation value 
of the electron velocity are high. In a close encounter of the 
antiproton on a molecular target the electron cloud might be 
more efficiently moved away from the projectile towards the 
other nuclei since there is in contrast to atoms always one pos- 
itive particle which is not neutralized by the antiproton. 
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Figure 3: (Color online) Cross sections for (a) single ionization and 
(b) excitation of molecular hydrogen by antiproton collisions for dif- 
ferent molecular orientations. Black pluses, orientationally averaged; 
red circles, (i); blue squares, (ii); green triangles, (Hi). The inset 
shows a sketch of the three orientations in the molecule-fixed frame. 



At high energies the single-ionization cross section for he- 
lium and 4 times the molecular hydrogen ion are both similar 
to the curve of the hydrogen atom while the results for molec- 
ular hydrogen are in good agreement with twice the curve for 
the hydrogen atom [cf. Fig.[Tfa)]. For these energies distant 
encounters are dominating the ionization process. Accord- 
ingly, details of the targets like the exact distribution of the 
positive charges become less important and the cross sections 
are mostly determined by the ionization potential. 

In Fig. [3] the dependence on the orientation of the molecu- 
lar axis with respect to the antiproton trajectory is presented. 
The cross sections er(0, $) for (a) ionization and (b) excita- 
tion are given for the three orthogonal orientations (0, $) = 
(i) (0, 0), (ii) (tt/2, 0), and (Hi) (tt/2, tt/2) (cf. the sketch in 
Fig. [3] and Ref. lfl7IO revealing the following results. First, 
the curves for the three different orientations generally dif- 
fer considerably especially around the maximum. Second, the 
calculation of only the parallel orientation (i) reproduces for 
energies above the maximum the orientationally-averaged re- 



sults with less than 3 % relative deviation. Third, for lower 
energies at which close collisions become dominant the con- 
sideration of the molecular geometry is inevitable. In contrast 
to the findings for the molecular hydrogen ion fl7il . the curves 
of orientation (i) are close to those of (ii) for ionization be- 
low the maximum and of (Hi) for excitation below 5 keV. In 
general, the differences among the cross sections for the three 
orientations are less pronounced than for the molecular hydro- 
gen ion. This might be due to the smaller internuclear distance 
and the two electrons of the hydrogen molecule making it a 
more spherical target. 

In conclusion, theoretical data are presented for single ion- 
ization and excitation of molecular hydrogen by antiproton 
impact for a wide energy range obtained with a two-electron 
time-dependent close-coupling method. The experimental 
data are in good agreement with the present calculations at 
high energies but are larger around the maximum. For en- 
ergies below the maximum the ionization cross section de- 
creases with decreasing energy much faster than in the cases 
for the hydrogen and helium atom but in a similar way as 
for the molecular hydrogen ion revealing the differences be- 
tween atoms and molecules. Furthermore, the importance of 
the molecular geometry and a full two-electron description is 
demonstrated. The present work should motivate new experi- 
mental efforts for molecular targets at low impact energies to 
confirm and further extend the gained insight. Additionally, 
it provides benchmark data for molecular collisions in general 
and for single ionization and excitation of molecular hydrogen 
by antiproton impact, in particular, which might be useful for 
the development of molecular modes. 

The method should be further exploited to extract also 
quantities like, e.g. electron-energy spectra or differential ex- 
citation cross sections. Such results may allow, especially at 
low energies, for the elimination of the diversity of results for 
the stopping power obtained in different experiments ifToL fl~3tl . 
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